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Abstract. Although quasi- two-dimensional organic superconductors such as n- 
(BEDT-TTF) 2 Cu(NCS)2 seem to be very clean systems, with apparent quasiparticle 
mean-free paths of several thousand A, the superconducting transition is intrinsically 
broad (e.g ~ 1 K wide for T c « 10 K). We propose that this is due to the extreme 
anisotropy of these materials, which greatly exacerbates the statistical effects of spatial 
variations in the potential experienced by the quasiparticles. Using a statistical 
model, we are able to account for the experimental observations. A parameter x, 
which characterises the spatial potential variations, may be derived from Shubnikov- 
de Haas oscillation experiments. Using this value, we are able to predict a transition 
width which is in good agreement with that observed in MHz penetration-depth 
measurements on the same sample. 
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The large number of Shubnikov-de Haas and de Haas-van Alphen oscillation 
experiments which have been carried out on crystalline organic superconductors 
demonstrates the high quality of these materials 0; oscillations are resolved down to 
~ 2 T [0, and apparent scattering times extracted by Dingle analysis |1], |3|] are well 
in excess of a picosecond, suggesting intralayer mean- free paths >1000 A §. Support 
for the cleanliness of the organics also comes from magneto-optical measurements of 
cyclotron- and Fermi-surface-traversal resonances, which again yield apparent scattering 
times >1 ps @, §. 

In spite of this, the superconducting transition at the critical temperature T c seems 
to be rather broad, whatever the measurement method used. Resistivity measurements 
are the most prone to complications || [?]], especially in an applied magnetic field. 
However, even when reliable perturbative techniques such as thermal conductivity [g] 
or GHz and MHz penetration-depth experiments || || |§ are employed in zero field, 
the transition has a significant width AT C . Figure [I] is a typical example; a k-(BEDT- 
TTF) 2 Cu(NCS)2 sample was placed in a coil forming part of a tank circuit oscillating 
at around 38 MHz (see Ref. |J). The superconducting to normal transition is observed 
because of the change from skin-depth to penetration-depth limited coupling of the 
sample to the MHz fields || |{|, which results in a shift in resonant frequency / of the 
tank circuit. For the purpose of making quantitative comparisons below, we choose two 
methods for defining the temperature (T) width of the transition. Firstly, by fitting the 
differential df/dT of the data to a Gaussian centered on T c uss = 9.38 K (Figure 
inset), a full width of AT c Gauss « 0.7 K is obtained. Alternatively, AT C may be defined 
using straight-line extrapolations (see Figure [I]), to give AT c lmcar m 0.9 K. Note that 
both methods yield AT C /T C ~ 0.1; it is likely that this significant intrinsic broadening 
of the transition region is responsible for the wide range of T c values quoted for k- 
(BEDT-TTF) 2 Cu(NCS) 2 [|, %. 

In this Letter, we show how a superconducting transition can be broadened due 
to charged impurities (or vacancies) randomly dispersed throughout a crystal, even 
when the potential associated with the impurity (or vacancies) is of extremely short 
range. When averaged over the coherence volume A v = ^03.^^02/6, where the £ 
are Pippard coherence lengths,!]] statistical variations in the density of impurities (or 
vacancies) lead to a spatially varying order parameter A and consequently a Gaussian- 
broadened transition. The effects of these statistical variations are shown to become 
more pronounced when the dimensionality of the superconductor is reduced, as is the 
case in quasi-two-dimensional crystalline organic superconductors such as k-(BEDT- 
TTF) 2 Cu(NCS) 2 . 

Our treatment is closely related to the Ginzburg criterion jR], [11], a quantitative 
guide to the circumstances under which mean-field theory can be expected to break 

|| As the following paragraph will show, we are dealing with a situation in which mean-field theory 
breaks down because of spatial variations of the order parameter. It is therefore inappropriate to use 
the Ginzburg-Landau temperature dependent coherence lengths to parameterise the coherence volume, 
as these diverge at T c . 
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Figure 1. MHz penetration data for a single crystal of k-(BEDT-TTF) 2 Cu(NCS) 2 , 
shown as resonant frequency / versus temperature T. The superconducting transition 
is the steeply sloping region between the more gentle variations characteristic of 
superconductivity (low T) and the normal state (high T); note that the complete 
transition region occupies a temperature range from around 8 K to 11 K. The dotted 
lines are extrapolations of the normal-state, transition-region and superconducting- 
state behaviour. The intersections of the extrapolations occur at 8.9 T and 9.8 T, giving 
r nncar _ 9 35 K ( mi d po int) and AT c lincar w 0.9 K. The inset shows the differential 
d//dT of the data (points) fitted to a Gaussian (curve) centered on y c Gauss — 9.38 K, 
with a full width of AT„ Gauss w 0.7 K. 



down due to fluctuations.^] The Ginzburg criterion has in the past been invoked to 



explain the broad superconducting transitions in the "high T c " cuprates |[12| . However 



in contrast to the situation in the cuprates, details of the bandstructures of the organic 
superconductors are often known to great precision |IJ, whereas the interpretation of 
heat capacity data (necessary to derive the Ginzburg criterion) in the latter systems is 
still somewhat contentious 0. We have therefore used an alternative statistical method 
to treat the spatial variations of the potential, based on a good quantitative knowledge 
of g(E F ), the density of quasiparticle states close to the Fermi energy E F . 

The introduction of impurities in an ideal metal leads to potential variations (of 
typical length scale R) and a finite scattering rate r _1 for normal quasiparticles jl], [T3|j . 
By contrast, the spatial extent £ of the superconducting wavefunction often greatly 
exceeds R, so that theoretical studies fail to yield a direct correlation between r _1 and 
% A useful introduction to the Ginzburg criterion is given in Section 5.1 of Reference uM. 
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the order parameter A ||1 5l.pl Nevertheless, nonmagnetic impurities can have an effect 
because the superconducting state is sensitive to changes in g(E F ). As an example of 
this, let us consider the weak-coupling BCS formula 

Ao-^oexpf--^], (1) 

where V is an interaction strength |15[]. If a small (local) fraction x of the host atoms or 
molecules are replaced by impurities, then the variation of g(E F ,x) can be written [|T^] 

g(E F ,x)=g(E F )+g'(E F )x. (2) 

If we assume that each impurity introduces one extra charge (— e), the derivative 
g'(E F ) = (dg(E F ,x)/dx) can be obtained from bandstructure calculations in the limit 
x -> 0. 

The number of sites available for impurity substitution within A v is n = A v /u v , 
where u v is the formula-unit volume. However, the local number of impurities within 
such a volume, m = xn, will be subject to statistical variations via the Binomial 
distribution (BD) |I 



, . x[l — x] n m n\ , . 

p{m,n) = — , 3 

ml [n — m\ ! 

where x is the mean of x. Below we consider relatively large values of n (43<n<10 10 ); the 
skewness factor r\ = 1/ ^Qx[l — x]n vanishes for large n, and the BD is well approximated 
by the normal distribution ||13||. Hence, the mean value of m becomes nx while its 



standard deviation is a(m) = Jnx[l — x]. In a "clean" metal x <C 1, leading to a 
standard deviation in x of 



y/x/n- ( 4 ) 

Using the fact that the standard deviation of a function of x is equal to the derivative 
of that function multiplied by the standard deviation of x, we insert Equations and f| 
into Equation |T|, yielding the standard deviation cr(A ) of A , 



g ( A o) _ 9'(E F ) xu v 
A ~ g{E F fV\ A v ■ 1 ] 



Applying the relation 2A = ak^Tc, where a ~ 3.52 in the weak-coupling BCS limit [15 
the broadening of the superconducting transition is 
AT C a(T c ) a(A ) 

T c T c A { } 

We now demonstrate the sensitivity of AT C to the dimensionality of the superconductor. 

Conventional three-dimensional (3D) superconductors have coherence lengths ~ 
10 2 - 10 4 A and u v ~ 20 A 3 0, yielding 10 4 <n<10 10 . In reasonably pure 
metals, we expect 10 _4 <a;<10~ 2 , while g(E F )V ~ 0.3; finally we use the fact that 
g'(E F )/g(E F ) — | for a free-electron model fl4l| . These figures yield transition widths 
10 _3 <AT C /T C <10 -7 , in reasonable agreement with observations |15|. These sharp 

+ Except in the case of magnetic impurities |fq] . 
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transitions are a consequence of the large size of the superconducting wavefunction, 
allowing inhomogeneities to be averaged out. 

In quasi-two-dimensional (Q2D) organic superconductors, the intralayer coherence 
lengths are ~ 10 2 A. However, the inter-layer coherence length £ 0z can be much less 
than the layer spacing a |7|]; in k-(BEDT-TTF) 2 Cu(NCS)4, £02 can be estimated to 
be ~ 0.3 A,Q whereas a m 16 A The consequence of this extreme anisotropy is 

that A v is replaced by a "coherence area" A a = vr£ x£oj//4, because superconducting 
wavefunctions do not extend out of the layers. In enumerating n, the unit cell must also 
be represented by an area w a (^ 106 A 2 in k-(BEDT-TTF) 2 Cu(NCS) 4 @). Assuming 
BCS values for the intralayer coherence lengths (£oz ~ (,oy ~ 76 A) jy|, we obtain 
n pd 43, much less than typical values of n in 3D systems. 

In the current context, k-(BEDT-TTF) 2 Cu(NCS) 2 has the considerable advantage 
that its intralayer bandstructure may be represented to good accuracy by the effective 
dimer model @, [19], p0[ 



E(k) = ±2 cos(^-)y / t 2 1 + t 2 2 + 2t cl t c2 cos(k c c) + 2t h cos(k h b) . (7) 

Here k\, and k c are the intralayer components of k and tb, ^ci and t c2 are interdimer 
transfer integrals |19|, EDI; the + and — signs result in the quasi-one- dimensional (Q1D) 
sheets and Q2D pocket of the Fermi surface respectively [I]]. Moreover, the Fermi- 
surface warping in the interlayer direction is rather small so that it may be ignored 
for the current purposes. Accurate de Haas-van Alphen and magnetic breakdown data 
constrain the parameters in Equation ^ rather tightly, leading to values t^ — 15.6 meV, 
t cl = 24.2 meV and t c2 = 20.3 meV §. 

Figure || shows N s , the areal quasiparticle density per layer, as a function of 
quasiparticle energy E for k-(BEDT-TTF) 2 Cu(NCS) 2 ; the curve has been derived using 
the parameters listed in the previous paragraph and Equation [7|. The position of the 
Fermi energy and the value of N s corresponding to the well-known magnetic breakdown 
(5 frequency Q, which encompasses the entire Fermi surface, are indicated by dotted 
lines. Whereas for many purposes the bands of K-(BEDT-TTF) 2 Cu(NCS)2are rather 
parabolic close to Ep, some curvature of A^ versus E is visible in Figure [|, resulting in 
a g{Ep) which varies with energy (Figure ||, inset). Using the results shown in Figure 
we obtain g'(E F )/g(E F ) ps 1.45. 

Whilst it is unlikely that k-(BEDT-TTF) 2 Cu(NCS) 2 is a weak-coupling BCS 
superconductor |], § |20], [21], p2| , Equation [I] is known to describe the functional 
dependence of T c on g(E F ) very well, as shown in the pressure-dependent experiments 
of Caulfield et al. 0. Using ^(^f)^ w 0.3 @ and 10- 4 <x<lQ- 2 we obtain 



0.015<AT C /T C <0.15, encompassing all known experimental data. 

In order to find a suitable value of x for the sample used in Figure [I], we turn to 
Reference |[L3|| , which shows that spatial variations in the potential experienced by the 

* The anisotropy is ~ 100 — 350 [|] |l^, |l^|, and the in-plane coherence lengths are estimated to be 
w 76 A 0,|18). 
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Figure 2. N s , the areal quasiparticle density per layer, as a function of quasiparticle 
energy E for k-(BEDT-TTF) 2 Cu(NCS) 2 ; the curve has been derived using the 
parameters listed in the text and Equation The position of the Fermi energy and the 
value of N s corresponding to the well-known magnetic breakdown f] frequency are 
indicated by dotted lines. The inset shows the variation of the quasiparticle density of 
states g(E), normalised to its value at E — E-p, as a function of quasiparticle energy. 



quasiparticles leads to broadening of Landau levels and hence damping of Shubnikov- 
de Haas and de Haas-van Alphen oscillations. This damping is parameterised by an 
effective Dingle temperature 



In the current context, F is the quantum oscillation frequency of the (3 breakdown orbit 
and F' = dF/dx; it is simple to show that F' = F pa 3920 T |]. The sample used in 
Figure [I] had T Da pa 0.42 K j| for the a-orbit Shubnikov-de Haas oscillations, which 
translates into a Dingle temperature of Trj ~ 0.76 K for the f3 magnetic breakdown 
orbit [|TJ. Substituting this value of T D into Equation || yields x ~ 0.0017, which can 
be used in Equation ^| to give AT C pa 0.6 K. This value is in good agreement with 
AT c Gauss pa 0.7 K extracted from the MHz experiments (see Figure [I|, inset) and close to 
the value AT c linear pa 0.9 K found using the alternative linear extrapolation method. 

One possible interpretation of the value of x ~ 0.0017 is that ~ 0.1 — 0.2% of the 
molecular sites in our crystals of k-(BEDT-TTF) 2 Cu(NCS)2 are in some way defective, 
functioning as "impurities" or "vacancies"; possible mechanisms might include neutral 
BEDT-TTF molecules, anions which are missing or in which the Cu ion possesses the 




(8) 
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wrong charge, or the incorporation of other molecular species from the growth process. 
The detection of such defects by other means is very difficult; all that can be noted at 
present is that such concentrations of defects are thought to be not unlikely [p3] . 

In summary, using a statistical model, we are able to account for the broadened 
superconducting to normal transitions observed in organic superconductors. Our 
model consistently explains the superconducting transition width and the Landau-level 
broadening in k-(BEDT-TTF) 2 Cu(NCS)4 using one parameter, x. Given the precise 
knowledge of the Fermi-surface topologies of many organic superconductors, and the 
availability of many samples of differing quality and effective dimensionality ||, our 
model may be useful in obtaining a more detailed understanding of the factors which 
influence the formation of the superconducting state. 
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